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The size of non-perturbative corrections to event shape observables is predicted to fall like 
powers of the inverse centre of mass energy. These power corrections are investigated for dif- 
ferent observables from e + e~-annihilation measured at LEP as well as previous experiments. 
The obtained corrections are compared to other approaches and theoretical predictions. Mea- 
surements of a s using power corrections are compared to conventional methods. 



1 Introduction 

The process of hadron production in e + e _ -annihilation is usually depicted by three phases. 
The first part called the perturbative phase is described by perturbative QCD calculations. The 
second phase called fragmentation or hadronisation is usually described using Monte Carlo based 
models. It is widely hoped that the influence of this phase on event shape observables can be 
described by analytical means: Power corrections. A third phase containing hadron decays is 
believed to be well under control. 

Power corrections arise from two different theoretical approaches: Renormalons and analyt- 
ical hadronisation models. This suggests a connection between the picture of a hadronisation 
phase and the theoretical idea of renormalons. 

The size of the correction due to the hadronisation process can be seen as a quality attribute 
for a specific observables. A small correction implies that this observable probes the parton 
structure more directly, allowing a more precise test of QCD predictions, e.g. measurements of 
a s . As power corrections use few parameters for describing the hadronisation phase, it is hoped 
that they lead to a better understanding of its influence. 



2 Simple Power Corrections 



2.1 Mean Values 

Means of infrared and collinear safe observables can be described by the sum of the perturbative 
part and the power correction term 

(/) = —/ ffda = (/ pert ) + </ pow ) (1) 

ftot J acr 

where the perturbative prediction in 2nd order can be written as 

with A and B being given numbersBi, fi being the renormalisation scale and &o = (33— 2Nf)/12n. 
To investigate the size and type of power corrections, 

(/pow) = p I" jp2~ (3) 
"cm "cm 

is used as a simple ansatz. It is useful to fix a s in these fits to get comparable power coefficients. 

The four mean event shapes in Figure |l] show qualitative agreement between the parton 
levels of the parton shower Monte Carlo and the second order part resulting from the fit. The 
corrections of all four means show \/E cm behaviour with only 1 — T having large x 2 due to 
inconsistent data. Although the Ci coefficients of (1 — T) and {B snm ) are not consistent with 
0, the dominant contributions comes from the 1/E cm term. (1 — T) and {B snm ) have larger 
corrections than (M^/E^ is ) and (B mSLX ). The numbers are given in Table [p. 

2.2 Cut Integrals and Higher Moments 

To investigate whether the non-perturbative correction of an event depends on specific values of 
an observable or not cut integrals (x) a b := J a b ^ dx were investigated. While (1 — T) Q 1 5 
shows dominantly 1/E cm correction like the full mean, (Af^/^ is ) 01 05 acquires a 1/E^ m be- 
haviour. (Table [l] and Figure ||a). 

Another way of emphasising different ranges is to investigate higher moments. The OPAL 
Coll. investigated the power dependence of the first three moments of 1 — T and the C- 
parameter @. 



(/pow) = a f \ (-jjr-) atn-ifjii) ~ at,(n) - ( b ■ In ^ + f + — J o? a (ji) 
n l \E cm J \ 2ir n J 
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was used as power term. a n _i being a non-perturbative parameter accounting for the contribu- 
tions to the event shape below an infrared matching scale fij, K = (67/18 — 7r 2 /6)Ca — 5iV#/9 
and (if = ACf/n. Beside a s these formulae contain a n _i as the only free parameters. The 
results shown for (1 — T) in Figure ^ and Table ^ indicate that the assumed power law of 1/E™ m 
for the n-th moment does work. 



Table 1: Coefficients from fitting Equations (pHa) to various observables with a s (Mz) = 0.118 fixed. 



Observable 


Ci [ GcV] 


C 2 [ GeV 2 ] 


x 2 

ndf 


Observable 


Ci [ GeV] 


C 2 [ GeV 2 ] 


x 2 

ndf 


(1-T) 
(K/^s) 

(^max) 


1.09 ±0.03 
0.67 ±0.12 
1.27 ±0.09 
0.33 ±0.06 


- 3.2 ±0.9 

- 1.9 ±0.7 
-10.3 ±4.4 

0.2 ±2.3 


72/38 
21/25 
21/19 
23/20 


(1 ~ ^/0.2...0.5 
\ ^vis / 0.1. ..0.5 


0.37 ±0.03 
0.08 ±0.03 


0.6 ±1.2 
9.2 ±0.9 


20/9 

8/5 
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Figure 1: Measured (1 — T), (M^/E^ is \, (B sum ) and {_B max )Baffl as a function of the centre of mass energy. On 
the lefthand-side the solid lines present the hadron level prediction of Jetset74, the dotted lines show the parton 
shower result. On the righthand-side the solid lines present the results of the fits with Eqs. (Q-j^), the dotted 
lines show the perturbative part only. The strong decrease of the {B sum ) prediction at low energies at the right 
is caused by the large and negative C2 coefficient which is poorly determined. 
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Figure 2: Results of fitting Eqs. (0,^,^) to Jetset predictions for mean values and higher moments of 1 — T. 

Table 2: Results of fitting Eqs. to Jetset predictions for mean values and higher moments of 1 — T. 



1 - T 


Mean 


2nd Moment 


3rd Moment 


a s (M z ) 


0.108 


0.1121 


0.1131 


«n-l 


0.380 


1.34 


5.7 


X 2 /ndf 


16 


4.8 


3.0 



3 Other Approaches 

To strengthen the need for power corrections in the description of data, a few other approaches 
were investigated: Fitting only the perturbative prediction including a 3rd order coefficient as 
free parameter yields very large x 2 an d a 3rd order coefficient C = 0(1000) for (1 — T) and 
(M^/E 2 is ). Thus a 3rd order calculation does not give an improved prediction on the energy 
dependence neglecting hadronisation effects. (Figure ||b). 

Leaving the renormalisation scale = [i 2 /E 2 m free, shows ambivalent results. While for 
(1 — T) the fit is reasonable and leads to a scale that is consistent with the one obtained 
from fits to distributions, the fit for (M^/E 2 is ) cannot describe the data. (Figure |||c). 

To obtain the functional type of power corrections needed, one can also fit with an arbitrary 
power p: 

(/pow) = -jJr- (5) 
For (1 — T) one gets p = 0.98 ± 0.1!)i in perfect agreement with the previous results. 

4 Determination of a s using Power Corrections 

The analytical power ansatz for non-perturbative corrections by Dokshitzer and Webber i'0 is 
used by DELPHI to determine a s from mean event shapes BQ. The power correction of this 
ansatz is given by 



it \ VI 
(Jpow) - a f ■ — 



ao(fJ-i) - a s (jx) - ( b ■ In ^ + ^ + 2b o ) a? s (fj ) 



(6) 



with [aj, K, af and &o as given above. In order to measure a s from individual high energy data 
the free parameter o.q has to be known. 

To infer uq a combined fit of a s and a$ to a large set of measurements at different energies! 
is performed. For E cm > Mz only DELPHI measurements are included in the fit. The resulting 
values of ao are summarised in Table ||. The extracted ao values are around 0.5 as expected int3. 
The numerical values are, however, incompatible with each other. So the assumed universality^ 
is not valid to the precision that is accessible from the data. Therefore «o is determined for 
(1 - T) and (Mfi/E^) individually The scale error is obtained from varying the renormalisation 
scale factor x M = iJ, 2 /E^ m from 0.25 to 4 and the infrared matching scale from 1 GeV to 3 GeV. 



After having fixed ao, the a s values corresponding to the high energy data points can 
be calculated from Eqs. (|I|j2|,^). a s is calculated for both observables individually and then 
combined with an unweighted average. The resulting a s values and the QCD expectation are 
shown in the leftmost plot of Figure |4j. 

The a s values follow the QCD expectation. Fitting a straight line to the energy dependence 
results in a slope, which agrees very well with the QCD expectation of a running a s between 
91 GeV and 183 GeV. In Figure [| and Table || this result is further-on compared to a s mea- 
surements obtained from distributions using Monte Carlo based hadronisation models. The 
comparison shows that the running measured from means using power correction and the one 
obtained from standard methods give consistent results and comparable errors. 

5 Prospects 

Many of the recent developments in the field of power corrections are not yet included in the 
presented experimental works. Beside the calculation of Milan factors, which will not influence 
the consistency of the experimental results much, some of the predicted power coefficients a/ 
were corrected in the last months til. First tests show that the consistency of measured <5o and 
a s values improve with these new predictions. 

Thus the experimental fits have to be repeated and with increasing theoretical understanding 
consistent experimental results hopefully will arise. 
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Figure 3: Fit results for cut variables (left), fitted 3rd order coefficient (middle) and fitted renormalisation 
scale (right). The continuous line represents the complete prediction. In the left plot the dashed line gives the 
perturbative contribution. In the middle plot it represents the 2nd order contribution to the full 3rd order curve. 



Table 3: Determination of ao- For E cln > Mz only DELPHI measurements are included in the fit. The first error 
is the statistical error from the fit, the second one is the scale error. 



Observable 


a 


a s (M z ) 


Aust MeV] 


X 2 /ndf 


(l-T) 
{M 2 h /El- ls ) 


0.531 ±0.012 ±0.003 
0.434 ±0.010 ±0.010 


0.1187 ±0.0017 ±0.0060 
0.1144 ±0.0012 ±0.0041 


235 ± 22 ± 70 
184 ± 13 ± 40 


42/23 
4.0/14 




Figure 4: Energy dependence of a 3 as obtained from mean event shapes (left) compared to a s obtained from 
distributions. The errors shown are statistical only. The band shows the QCD expectation of extrapolating the 

world average to other energies. 



Table 4: Running of a s (da/dE crn ) between 91 GeV and 183 GeV. 



Means plus 
Power corr. 


1 

Monte Carle 


Distributions pi 
) based hadroni 
NLLA 


us 

sation models 
0(a2)+NLLA 


QCD pred. 


-0.13 ±0.05 


-0.13 ±0.05 


-0.15 ±0.05 


-0.13 ±0.04 


-0.13 • 10~ 3 GeV" 1 



